Nonuniversal mound formation in nonequilibrium surface growth 
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We demonstrate, using well-established nonequilibrium limited-mobility solid-on-solid growth mod- 
els, that mound formation in the dynamical surface growth morphology does not necessarily imply 
the existence of a surface edge diffusion bias ("the Schwoebel barrier"). We find mounded mor- 
phologies in several nonequilibrium growth models which incorporate no Schwoebel barrier. Our 
numerical results indicate that mounded morphologies in nonequilibrium surface growth may arise 
from a number of distinct physical mechanisms, with the Schwoebel instability being one of them. 
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In vacuum deposition growth of thin films or epitax- 
ial layers (e.g. MBE) it is common 111 to find mound 
formation in the evolving dynamical surface growth mor- 
phology. Although the details of the mounded morphol- 
ogy could differ considerably depending on the systems 
and growth conditions, the basic mounding phenomenon 
in surface growth has been reported in a large num- 
ber of recent experimental publications [jj]. The typi- 
cal experiment [lj monitors vacuum deposition growth 
on substrates using STM and/or AFM spectroscopies. 
Growth mounds are observed under typical MBE-type 
growth conditions, and the resultant mounded morphol- 
ogy is statistically analyzed by studying the dynamical 
surface height h(r, t) as a function of the position r on 
the surface and growth time t. Much attention has fo- 
cused on this ubiquitous phenomenon of mounding and 
the associated pattern formation during nonequilibrium 
surface growth for reasons of possible technological inter- 
est (e.g. the possibility of producing controlled nanoscale 
thin film or interface patterns) and fundamental interest 
(e.g. understanding nonequilibrium growth and pattern 
formation) . 

The theoretical interpretation of the mounding phe- 
nomenon has often been based jjj on the step-edge diffu- 
sion bias H or the so-called Schwoebel barrier || effect 
(also known as the Ehrlich-Schwoebel ||, or ES, barrier). 
The basic idea of the ES barrier-induced mounding (often 
referred to as an instability) is simple : The ES effect pro- 
duces an additional energy barrier for diffusing adatoms 
on terraces from coming "down" toward the substrate, 
thus probablistically inhibiting attachment of atoms to 
lower or down-steps and enhancing their attachment to 
upper or up-steps; the result is therefore mound forma- 
tion because deposited atoms cannot come down from 
upper to lower terraces and so three-dimensional mounds 
or pyramids result as atoms are deposited on the top of 
already existing terraces. 

The physical picture underlying mounded growth un- 
der an ES barrier is manifestly obvious, and clearly the 
existence of an ES barrier is a sufficient condition || 
for mound formation in nonequilibrium surface growth. 



Our interest in this paper is to discuss the necessary con- 
dition for mound formation in nonequilibrium surface 
growth morphology — more precisely, we want to ask 
the inverse question, namely, whether the observation of 
mound formation requires the existence of an ES bar- 
rier. Through concrete examples we demonstrate that 
the mound formation in nonequilibrium surface growth 
morphology does not necessarily imply the existence of an 
ES barrier, and we contend that the recent experimental 
observations of mound formation in nonequilibrium sur- 
face growth morphology should not be taken as defini- 
tive evidence in favor of an ES barrier-induced univer- 
sal mechanism for pattern formation in surface growth. 
Mound formation in nonequilibrium surface growth is a 
non-universal phenomenon, and could have very differ- 
ent underlying causes in different systems and situations, 
with the Schwoebel instability being one particular mech- 
anism (among many) for the mounded morphology. 

Before presenting our results we point out that the 
possible nonuniversality in surface growth mound forma- 
tion (i.e. mounds do not necessarily imply an ES barrier) 
has recently been mentioned in at least two experimen- 
tal publications |4||| where it was emphasized that the 
mounded patterns seen on Si |H and GaAs §, InP § 
surfaces during MBE growth were not consistent with 
the phenomenology of a Schwoebel instability. In two 
other recent experimental publications || mound forma- 
tion during semiconductor surface growth (Ge, GaAs) 
was carefully analyzed using the prevailing Schwoebel 
instability phenomenology with a conclusion not very 
dissimilar from that in ref. In particular, the ES 
barrier-based analyses of the experimental data in both 
the papers in refs. || produced rather weak Schwoebel 
effects in both experiments, leading to the conclusion in 
both experiments that the Schwoebel instability in all 
likelihood is playing a small to negligible role in the ob- 
served mound formations in refs. |6). Very recently, ex- 
perimental observations of striking mound formation JtJ 
in Au and MgO vapor deposition growth have been in- 
terpreted without invoking any ES barrier effect. Thus, 
the observed mound formation in the nonequilibrium sur- 
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face growth in refs. ||-|7[ is interpreted essentially with- 
out invoking any key role being played by an ES barrier 
whereas the mounded growth morphologies in refs. 
have mostly been interpreted as arising essentially due to 
a Schwoebel instability. Thus the inevitable conclusion 
from recent experimental observations l^-Q is that the 
mound formation in surface growth does not necessarily 
arise from the universal mechanism of a Schwoebel insta- 
bility §, but may be caused by different non-universal 
mechanisms in different experimental situations. The 
purpose of the current article is to explore this nonuniver- 
sality in the mound formation in some detail using sim- 
ple solid-on-solid (SOS) growth models where the kinetic 
mechanisms leading to the mounded morphologies are 
explicitly obvious, and therefore compelling conclusions 
can be drawn about the precise physical mechanism pro- 
ducing the mounds. A direct comparison between exper- 
imental results and our rather simplistic limited mobil- 
ity nonequilibrium SOS models, however, is unwarranted 
due to the extreme simplicity in the growth and diffusion 
rules in our models — our models do serve the purpose 
of explicitly demonstrating the fact that mounded mor- 
phologies can arise without any ES barriers whatsoever. 

There have been two proposed mechanisms in the lit- 
erature which lead to mounding without any explicit ES 
barrier: One of them invokes Q a preferential attachment 
to up-steps compared with down-steps (the so-called 
"step-adatom" attraction), which, in effect, is equivalent 
to having an ES barrier because the attachment prob- 
ability to down-steps is lower than that to up-steps ex- 
actly as it is in the regular ES barrier case |2|,[| — we 
therefore do not distinguish it from the ES barrier mech- 
anism, and in fact, within the simple growth models we 
study, these two energetic mechanisms are physically and 
mathematically indistinguishable. The second mound- 
ing alternative P], which is a purely topologic-kinetic 
effect, is the so-called edge diffusion induced mounding, 
where diffusion of adatoms around cluster edges is shown 
to lead to mound formation during nonequilibrium sur- 
face growth even in the absence of any finite ES bar- 
rier. One of the concrete examples we discuss below, 
the spectacular pyramidal pattern formation (Fig. 3(c)) 
in the 2+1 dimensional (d) noise reduced Wolf- Villain 
(WV) model |fiof , arises from such a nonequilibrium edge 
diffusion effect (perhaps in a somewhat unexpected con- 
text). We also demonstrate, using the WV model and the 
Das Sarma-Tamborenea (DT) model [[ll|, that mound 
formation during nonequilibrium surface growth is, in 
fact, almost a generic feature of limited mobility solid- 
on-solid discrete growth models JlO|-p^|, which typically 
have comparatively large values of the roughness expo- 
nent (a) characterizing the growth morphology. We 
find that a large roughness exponent coupled with atom- 
istic solid-on-solid growth almost invariably leads to vi- 
sually mounded growth morphology. Below we demon- 
strate that mound formation in surface morphology aris- 



ing from this generic "large a" effect (without any ex- 
plicit ES barrier) is often qualitatively virtually indis- 
tinguishable from that in growth under an ES barrier. 
Mound formation in the presence of strong edge diffu- 
sion H (as in the d=2+l WV model in Fig. 3) is, on the 
other hand, morphologically quite distinct from the ES 
barrier- or the large a- induced mound formation. 

Our results are based on the extensively studied |l^] 
limited mobility SOS nonequilibrium WV and DT 
p"l| growth models. Both models have been widely stud- 
ied p"2[ in the context of kinetic surface roughening in 
nonequilibrium solid-on-solid epitaxial growth — the in- 
terest in and the importance of these models lie in the fact 
that these were the first concrete physically motivated 
growth models falling outside the well-known Edwards- 
Wilkinson-Kardar-Parisi-Zhang [[t2| generic universality 
class in kinetic surface roughening. Both models involve 
random deposition of atoms on a square lattice singular 
substrate (with a growth rate of 1 layer/sec. where the 
growth rate defines the unit of time) under the SOS con- 
straint with no evaporation or desorption. An incident 
atom can diffuse instantaneously before incorporation if 
it satisfies certain diffusion rules which differ slightly in 
the two models. In the WV model the incident atom can 
diffuse within a diffusion length I (which is taken to be 
one with the lattice constant being chosen as the length 
unit, i.e. only nearest-neighbor diffusion, in all the results 
shown in this paper — larger values of I do not change our 
conclusions) in order to maximize its local coordination 
number or equivalently the number of nearest neighbor 
bonds it forms with other atoms (if there are several pos- 
sible final sites satisfying the maximum coordination con- 
dition equivalently then the incident atom chooses one of 
those sites with equal random probability and if no other 
site increases the local coordination compared with the 
incident site then the atom stays at the incident site). 
The DT model is similar to the WV model except for 
two crucial differences: (1) only incident atoms with no 
lateral bonds (i.e. with the local coordination number 
of one — a nearest-neighbor bond to the atom below is 
necessary to satisfy the SOS constraint) are allowed to 
diffuse (all other deposited atoms, with one or more lat- 
eral bonds, are incorporated into the growing film at their 
incident sites); (2) the incident atoms move only to in- 
crease their local coordination number (and not to maxi- 
mize it as in the WV model) — all possible incorporation 
sites with finite lateral local coordination numbers are 
accepted with random equal probability. Although these 
two differences between the DT and the WV model have 
turned out to be crucial in distinguishing their asymp- 
totic universality class, the two models exhibit very sim- 
ilar growth behavior for a long transient pre-asymptotic 
regime. It is easy to incorporate |l3] | an ES barrier in the 
DT (or WV) model by introducing differential probabili- 
ties P u and Pi for adatom attachment to an upper and a 
lower step respectively — the original DT model pT| has 
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P u = Pi, and an ES barrier can be explicitly incorporated 
Jl3| in the model by having Pi < P u < 1 . We call this sit- 
uation HI the DT-ES model (we use P u = 1 throughout 
with no loss of generality). We also note, as mentioned 
above, that within the DT-ES model the ES barrier [gj| 
(Pi < P u ) and the step-adatom attraction fj (P u > Pi) 
are manifestly equivalent, and we therefore do not con- 
sider them as separate mechanisms. We note also that 
in some of our simulations below we have used the noise 
reduction technique p^p^ ] which have earlier been suc- 
cessful in limited mobility growth models in reducing the 
strong stochastic noise effect through an effective coarse- 
graining procedure. All three models described above are 
studied in both one-dimensional substrate (d=l+l) and 
two-dimensional substrate (d=2+l) systems with peri- 
odic boundary conditions being used in all simulations. 
Detailed descriptions of DT and WV models are available 
in the literature |lCj-0. 

In Fig. 1 and 2 we present our d=l+l growth simu- 
lations, which demonstrate the point we want to make 
in this paper. We show in Fig. 1 the simulated growth 
morphologies at three different times for four different 
situations, two of which (Fig. 1(a), (b)) have finite ES 
barriers and the other two (Fig. 1(c), (d)) do not. 
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FIG. 1. Dynamical morphologies at 10 2 , 10 4 and 10 6 
monolayers (ML) for (a) DT-ES with Pi = 0.5, P u = 1; (b) 
DT-ES with Pi = 0.9, P u = 1; (c) DT; and (d) WV models. 

The important point we wish to emphasize is that, while 
the four morphologies and their dynamical evolutions 
shown in Fig. 1 are quite distinct in their details, they 
all share one crucial common feature: they all indicate 
mound formation although the details of the mounded 
morphologies and the controlling length scales are obvi- 
ously quite different in the different cases. Just the mere 
observation of mounded morphology, which is present in 



Figs. 1(c), (d), thus does not necessarily imply the exis- 
tence of an ES barrier. To further quantify the mounding 
apparent in the simulated morphologies of Fig. 1 we show 
in Fig. 2 the calculated height-height correlation func- 
tion, H(r) ~ (h(x.)h(r + x)) x , along the surface for two 
different times. 
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FIG. 2. The height-height correlation function H(r) at 
10 2 ML (main plots) and 10 4 ML (insets) corresponding re- 
spectively to the morphologies in Fig. 1. 

All the calculated H(r) show noisy oscillations as a func- 
tion of r, which implies mound formation (corresponding 
to the noisy mounded morphologies of Fig. 1). It is 
indeed true that the presence of considerable stochastic 
noise associated with the deposition process in the DT, 
WV models make the H (r)-oscillations quite noisy, but 
the important feature to note in Fig. 2 is that the qual- 
itative oscillatory nature of H(r) in situations with or 
without an ES barrier is essentially the same. Thus, 
the mound formation, although noisy, is qualitatively 
similar with or without an ES barrier in Figs. 1 and 
2. We have explicitly verified that such growth mounds 
(or equivalently H(r) oscillations) are completely ab- 
sent in the growth models |0| which correspond to the 
generic Edwards- Wilkinson-Kardar-Parisi-Zhang univer- 
sality class, and arise only in the DT, WV limited mobil- 
ity growth models which are known to have large value 
of the roughness exponent a arising from (linear or non- 
linear) surface diffusion processes. In fact, the effective 
a in the DT, WV models is essentially unity (in 

d=l+l), which is the same as what one expects in a naive 
theoretical description of growth under the ES barrier 
(although the underlying growth mechanisms are com- 
pletely different in the two situations). We believe that 
any surface growth involving a "large" roughness expo- 
nent (0.5<a^l) will invariably show "mounded" mor- 
phology independent of whether there is an ES barrier in 
the system or not. We contend that this effectively large 
a is the physical origin for the mounded morphology in 
semiconductor MBE growth where one expects the sur- 
face diffusion driven linear or nonlinear conserved fourth 
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order (in contrast to the generic second order) dynami- 
cal growth universality [O] class to apply which has the 
asymptotic exponent : a (d=l+l) * 1; a (d=2+l) w 
0.67 (nonlinear), 1 (linear). One recent experimental pa- 
per ||, which reports the observation of mounded GaAs 
and InP growth with a w 0.5 — 0.6, has explicitly made 
this case, and other recent reported mound formations || 
in semiconductor MBE growth are also consistent with 
our contention that mounds may arise from a large ef- 
fective roughness exponent rather than a Schwoebel in- 
stability. Two very recent experimental publications |7]] 
have reached the same conclusion in non-semiconductor 
MBE growth studies — in these recent publications (?]] 
spectacular mounded surface growth morphologies have 
been interpreted on the basis of the fourth order con- 
served growth equations [^0[-p|. The crucial message 
of our simulated d=l+l growth morphologies in Fig. 1 
and 2 is the fact that the mound formation with [Figs. 
1 (a),(b) and 2 (a),(b)] and without ES barrier [Figs. 1 
(c),(d) and 2 (c),(d)] are qualitatively similar, and there- 
fore the mere observation of a mounded morphology does 
not necessarily imply a Schwoebel instability. 




FIG. 3. Morphologies from the (a) DT-ES with P t = 0.5, 
P u = 1; (b) DT; and (c) noise reduced WV models. 

Finally, in Figs. 3 and 4 we present our results for 
the physically more relevant d=2+l noncquilibrium sur- 
face growth. In Fig. 3(a)-(c) we show the growth mor- 
phologies for the DT-ES, DT, and the noise-reduced WV 



model, respectively whereas in Fig. 4 we show the scaled 
height-height correlation function for the mounded mor- 
phologies depicted in Fig. 3. It is apparent that all 
three models (one with an ES barrier and the other two 
without) have qualitatively similar oscillations in H(r) 
indicating mounded growth, and the differences in the 
mounding between the growth models are purely quanti- 
tative. Again, the important point is that mounded mor- 
phologies with and without ES barriers manifest similar 
oscillating in H(r), indicating that such an oscillatory 
height-height correlation function by itself does not es- 
tablish a Schwoebel instability. Thus we come to the 
same conclusion: mound formation, by itself, does not 
imply the existence of an ES barrier; the details of the 
morphology obviously will depend on the existence (or 
not) of an ES barrier. 
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FIG. 4. The scaled H(r) correlation functions correspond- 
ing to the morphologies shown in Fig. 3. (ro = mound ra- 
dius.) 

We note that the effective values of the roughness expo- 
nent are very similar in Fig. 3(a) and (b) i.e. with and 
without an ES barrier, both being approximately a ~ 0.5 
(far below the asymptotic value a « 1 expected in the ES 
barrier growth — we have verified that this asymptotic 
a 1 is achieved in our simulations at an astronomically 
long time of 10 9 layers). The most astonishing result we 
show in Fig. 3 is the spectacular pyramidal mound for- 
mation in the d=2+l noise reduced WV model (without 
any ES barrier), which has not earlier been reported in 
the literature. The strikingly regular pyramidal pattern 
formation (Fig. 3(c)) in our noise reduced WV model in 
fact has a magic slope and strong coarsening behavior. 
The pattern is very reminiscent of the theoretical growth 
model studied earlier in ref. 15 in the context of nonequi- 
librium growth under an ES barrier where very similar 
patterns with slope selection were proposed as a generic 
scenario for growth under a Schwoebel instability. In our 
case of the noise reduced d=2+l WV model of Fig. 3(c), 
there is no ES barrier, but there is strong cluster-edge 
diffusion. This strong edge diffusion (which obviously 
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cannot happen in 1+1 dimensional growth) arises in the 
WV model (but not in the DT model) from the hop- 
ping of adatoms which have finite lateral nearest neigh- 
bor bonds (and are therefore the edge atoms in a clus- 
ter). This edge diffusion (discussed in entirely different 
contexts in || ) leads to an "uphill" surface current in the 
111 direction, which leads to the formation of the slope- 
selected pyramidal patterned growth morphology. While 
noise reduction enhances the edge current strengthening 
the pattern formation (the uphill current is extremely 
weak in the ordinary WV model due to the strong sup- 
pression by the deposition shot noise), our results of Fig. 
3 estabish compellingly that the WV model in d=2+l 
is, in fact, unstable (uphill current) in contrast to the 
situation in d=l+l. Thus, the WV model belongs to 
totally different universality classes in d=l+l and 2+1 
dimensions! We mention that in (unphysical) higher (e.g. 
d=3+l, 4+1, etc.) dimensions, the WV model would be 
even more unstable, forming even stronger mounds since 
the edge diffusion effects will increase substantially in 
higher dimensions due to the possibility of many more 
configurations of nearest-neighbor bonding. Such unsta- 
ble growth in high dimensional (d > 2+1) WV model 
has earlier been reported [jl6| in the literature without 
any physical explanation. We have therefore provided 
the explanation for the long-standing puzzle of an insta- 
bility in high-dimensional (d > 2+1) WV model simu- 
lations which were reported |l6| in the literature some 
years ago. More details on this phenomenon will be pub- 
lished elsewhere |l7| ]. 

In conclusion, we have shown through concrete ex- 
amples that, while a Schwoebel instability is certainly 
sufficient to cause mounded surface growth morphol- 
ogy, the reverse is not true : an ES barrier is by no 
means necessary to produce mounds, and mound forma- 
tion in nonequilibrium surface growth morphology does 
not necessarily imply the existence of a Schwoebel in- 
stability. In particular, we show that a large roughness 
exponent (without any ES barrier) as in the fourth or- 
der conserved growth universality class [p~o|— produces 
mounded growth morphologies which are indistinguish- 
able from the ES barrier effect. Any experimentally ob- 
served mounded morphology therefore requires a careful 
and detailed quantitative analyses to determine the 
physical mechanism (e.g. ES barrier, edge diffusion, large 
roughness exponent without any ES barrier) underlying 
its cause — in particular, the existence of a mounded 
growth morphology by itself may not imply the existence 
of any significant Schwoebel barrier. 
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